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Abstract 

Using the modified spin-wave theory we calculate static and dynamic spin 
structure factors in spin-liquid state of the J1-J2 model. The spectrum of 
excitations in the vector channel is discussed. The developed technique can 
also be applied to the t — J model describing copper oxide superconductors. 
PACS numbers: 75.50.Ee, 75.10.Jm, 75.10.-b 



1 



Typeset using REVTgX 



Magnetically disordered states in quantum spin models are of considerable interest. Much 
of this interest stems from the connection of this problem to high-T c superconductivity. The 
ground state of an undoped compound has long range antiferromagnetic order, which is well 
described by the Heisenberg model and has been studied by numerous methods.! However 
introducing a small number of holes leads to destruction of long range order. 

Destruction of long range order can be studied by introducing some frustration into the 
Heisenberg model. This simplified mechanism of the destruction is certainly far from the 
reality for cuprates, but the investigation of this relatively simple model allows us to develop 
approaches to more realistic models. From this point of view the main purpose of the present 
work is to develop a technique applicable to the t — J model. 

We will focus on the simplest frustrated Heisenberg model. It is the J\-Ji model defined 

by 

H = J 1 J2S r -S r , + J 2 J2 S r -S r ,. (1) 

NN NNN 

In this Hamiltonian, the J\ term describes the usual Heisenberg interaction of the nearest 
neighbor spins (S = |) on a square lattice, while the J 2 term introduces a frustrating 
interaction between the next nearest neighbor sites. For convenience, we set J\ — 1 and 
denote a = J%j J\ 

For small a, the ground state is Neel ordered. For large a the system is decomposed 
into two Neel ordered sublattices which, however, have the same quantization axis. This 
is so called collinear state. Whether or not the Neel and collinear states are separated in 
parameter space by a spin liquid state has been a subject of many discussions. Besides 
the spin-wave calculationsi~i, the model has been studied by the Schwinger boson mean 
field theory!, analysis of small lattices!, a series expansion!, a mean field theory of bond 
operators^, and other methods0. The majority of these works including most recent exact 
finite lattice dioganalizations@ and series expansionO indicate the existence of a different 
intermediate phase for 0.4 < a < 0.6. Therefore in this work we assume that the system 
undergoes a second order quantum phase transition at a = a c ~ 0.4 from the Neel to a 
spin-liquid state. We are interested in this spin liquid phase close to the transition point, 
so a > a c and a — a c -C 1. In our previous workB for description of this state we used 
modified spin-wave theory originally suggested by Takahashilil for the Heisenberg model at 
nonzero temperature. In the present work we use same approach. 
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We remind the reader that for spin-wave description it is convenient to use Dyson- Maleev 
transformational (see also review0) for localized spin S = 1/2, 

Sf = a], = (2S — a\ai)a h 

Si = S — ajai, for I G up sublattice; (2) 
£*m = ^mi S m = b m (2S — 6 m 6 m ), 

= —S + b^ m b m , for m e down sublattice, 

and Fourier representation for a\ and b m 

^ = J|E^ kri « k (3) 

' k 

E e ' krm& k. 

Here N is number of sites on square lattice. The summation over k, here and everywhere 
below is restricted to the inside magnetic Brillouine zone (\k x \ + \k y \ < n). There are two 
simple ways to find an effective Hamiltonian quadratic in the operators a and b . First way is 
just dropping the quartic terms, and this is the linear spin- wave theory (LSWT). Second way 
corresponds to the mean field treating of the quartic terms a^ab^b — > (a)a)b'b+ (crb')ab+ 
and this is mean field spin- wave theory (MFSWT). Both approximations give very similar 
results (see e. g. Ref.@), Therefore as soon as we believe that spin-liquid phase exists, it 
does not matter which approximation is used for calculations of the properties of this phase. 
However LSWT gives the value of a c ~ 0.4 very close to that found from exact numerical 
computations. Therefore in the present work we will use linear spin-wave theory. 

Using Bogoliubov transformation one can easily find spectrum of spin-waves and stag- 
gered magnetization of sublattice in the Neel state 



^ k = 2V[l-«(l-^ k )] 2 -7^ (4) 

- i/c \ i 1 2 1 
m = \(S Z )\ = 1 --^ 



As usually we have defined 



7k = - (cos k x + cos k y ) and 7/k = cos k x cos k y . 



At a = the staggered magnetization m = 0.3. It decreases with increasing of a and 
vanishes at a = a c ~ 0.4 where the system undergoes a second order transition into the 
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liquid state. Zero sublattice magnetization means that the ground state is a condensate of 
many spin waves a k and 6k- To describe the emerging phase we must take into account 
their non-linear interaction. We cannot do this exactly. However, there is an approximate 
methodEi originally suggested by Takahashilil for the Heisenberg model at nonzero tempera- 
ture. Following Takahashi, we impose an additional condition that sublattice magnetization 
is zero 

(S z u ~S z d ) = {\- a[a, + | - blb m ) = 0, (5) 

where u and d are the spin up and down sublattices. The constraint (|5|) gives an effective 
cutoff of unphysical states in Dyson-Maleev transformation. This question is discussed in 
the paperQ. 

The constraint @ is introduced into the Hamiltonian via a Lagrange multiplier |i/ 2 . 
Now we must diagonalize 

H v = H LSW T-\v 2 {S z u -Sl)^ (6) 

k 

where A k = 1 — ct(l — ^ k ) + |z/ 2 . The simple (linear) second term in (||), taken together 
with Eq. (d), takes account of non-linear interaction of spin waves. Diagonalizing Eq. (||) 
by Bogoliubov transformation 

a k = U k a k + Vk/3l k , (7) 
6-k = Vk4 + f/ k /3_ k , 

we get the spectrum of excitations 



= 2JAI - 7k \ (8) 



This spectrum has a gap u^j 1 + y| ~ z/, so the meaning of Lagrange multiplier is elucidated. 
Taking also into account that in thermal equilibrium 

n k = (a k a k ) = = 1 (9) 

exp(^ k /T) - 1 



we get from (H) the equation for v 



2 ^ A k 



= 1"T7E— (l + 2n k ). (10) 

k i^Vk 
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The spin-wave velocity does not vanish at critical point c = y2(l — 2a) ~ 0.7. Therefore 



for v, k <C 1 the spectrum is of the form u u ^ = \/v 2 + c 2 /c 2 and equation ([T0| ) can be re- 
writtenQ as 



m H > 



1 1 , (^/v 2 + c 2 k r 

—, / = c °th — 

ck ^v 2 + c 2 k 2 \ 2T 



0, (11) 



where m < is given by Eq.(^j). We would like to stress that that the condition c / 
is not crucial for the validity of the method. Moreover for t — J model, which we are 
mainly interested in, the speed vanishes, or even c 2 < 0. Nevertheless the method wor 



For J\ — J2 model c > 0, and we will use the simplification — > \Jv 2 + c 2 A; 2 , but from 
comparison with exact numerical solution of equation fllOD we know that it is valid only 
for very small v and T: u, T i. Solution of equation (|TTD at zero temperature is 
straightforward 

Vq = Tc\m\c 2 = nBc 2 (a — a c ) ~ 3.7(a — a c ). (12) 

We took into account that near critical point m = B(a c — a) with slope B w 2.4 according to 
LSWT and most recent exact finite lattice dioganalizationsEl and series expansion0. Index 



in vq indicates that it is a gap at zero temperature. Equation (11) can be also written a; 



IT In [2 sinh ^ J = -nmc 2 = v . (13) 

Similar equation has been obtained in the Ref.El for the non-linear a-model in the limit 
when the number of components of the order parameter M = 00. Solution of equation (|i~3"D 
at low and high temperature looks like 

v^v (l + e~ uo/T ), T<z/ , (14) 
v ^QT + — p=f , T > v , 

where 6 = 2 In = 0.9624. We follow the notations of Ref R 

Now we can proceed to the calculation of spin structure factor in spin-liquid state. It is 
quite similar to the Takahashi's calculation for Heisenberg model at nonzero temperature^. 
Using Dyson- Maleev representation @, constraint @ ((aj(t)ai(t)) = (6] n (t)fe m (t)) = |), and 
mean field procedure for averaging of the quartic terms ((a^ab^b) — > (a)a)(tfb) + (a/b<)(ab)) 
one finds 
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(s m (t)s,(o)} = -{bi^aimibu^m, 

(S P (f)S,(0)) = (4(t)aK0))(aK^)af(0)>, 



(15) 



where Z, /' G sublattice, and m G down sublattice. Further calculation gives 

(Si(t)Sj(0)) = f 2 (t,Tij) — -5 2 (t,Tij), i,j£same sublattice, 
(Si(t)Sj(0)) = — g (t, Tij), i,j G different sublattices. 



(16) 



where 



2 ^ ikr 



E* 



e~^ k '(l + ri k ) + e w n k 



r ) = ;| £ e<kr [e"^ k '(l + 



n k ) - e*"" kt nk 



<?(*,r) 



2 e ^kr 7k 



1 + n k ) + e^ k 'n k 



(17) 



^From (|16|) and ([H]) one can easily find explicitly spin-spin correlator at large separation 



(m, n), r = \A« 2 + 



ir 



( 2u \ 1/iiTcr) 2 , 1 < r < cz//T 2 ; 

S t,r)S(t,0)) = (-ir + "exp( r) x A j ' 7 

V c / 1^ 2T/(7rc 3 z/r), r > cu/T 2 . 

So the magnetic correlation length equals £m = 0.5c/ v. 

Using (|16|) and (|T7| ) we can also find static and dynamic spin structure factors 



S M (q) = 5>^ r (S(t,r)-S(t,0)> 



N 



£ 



1 . AcAc+q - 7k7k+q . Q w-. . n 

t H -(1 + 2n k )(l + 2n k+q ) 



^k^i/k+q 



(18) 



(19) 



S M (w,q) = / — e'"'i:e">-'(S(i,r).S(0,0)) 



-E 

k 



1 AA+q ~ 7k7k+q 
"4 



(20) 



k-W^Vk+q 

x [(1 + n k )(l + n k+q )<5(^ - u vk - ^ k+q ) + n k n k+q 5(u; + uj vk + c^, k+q )] + 



+ 



1 ^A+q - 7k7k+q 
4 ^z/k^i/k+q 



2n k (l + n k+q )<!>(a; + u uk - ^ k+q ) ^ . 



We would like to stress that the dynamic structure factor Sm(u;, q) at T = contains 
only two-quasiparticle intermediate states. It is similar to the situation in one dimensional 
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Heisenberg chain where elementary excitation is spinon. We will discuss this point in con- 
clusion. Now let us calculate the structure factors SW(q) and Sm(^,<1) for q = and 
q = Q = (±7r, ±7r). Taking into account Eq.flD we find from (^0|) 

S „(M) = iMk*d+*) = ^H»<(" ,T T U/n T<< "°' (2D 

N k 7TC 2 L04T 2 1 > T > Z/ . 



One can easily obtain static factor Sm(0) integrating ( 2l|) by u. The calculation at q = Q 
is also very simple. It gives 

S M (Q)«-^lni, 

7TC Z Z/ 



with function F defined as 



v 



dx I 4T/t/ exp(-z/ /T), T < z/ ; 



/t x sinh 2 (:r/2) ) 1.67, 1 > T > u . 



The (^-distribution at q = and q = Q contains infinitely sharp 5(o>)-function and step 
function 6(±uj — 2u). However at q ^ 0, Q the cu-distribution is smooth. 
Average energy E = (H) can be easily calculated using Eqs. (|l|) and fllj 



E/iV = 2[-( 7 2 (0,r 1 ) + «/ 2 (0,r 2 )], (24) 

where ri = (1, 0) and r 2 = (1, 1). As an example, at Fig.l, we present the plot of energy as 
a function of temperature at a ~ a c ~ 0.4. The value of a — a c is chosen in such a way that 
the gap at zero temperature equals vq = 0.05. In this case our calculation gives for energy 
per site at zero temperature the value E(T = 0)/N = —0.513. This is slightly above the 
value E/N ~ -0.520 found in Ref.i by finite lattice dioganalization and in Ref.li3 by series 
expansion. 



Discussion 

In the present work using modified spin-wave theory we calculated static and dynamic spin 
structure factors in spin-liquid state of J\-Ji model. We believe that similar technique can 
be applied to the t — J model describing copper oxide superconductors. Let us point out 
the strong and weak sides of the modified spin-wave theory. It gives a very simple descrip- 
tion of the spin liquid state, and this is definitely a strong side. However the description 
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explicitly violates the rotational symmetry, and this is drawback of the approach. We would 
like to stress that violation of exact symmetry quite often appears in approximate descrip- 
tion of a strongly interacting system. For example Hartree-Fock method in atoms violates 
gauge-invariance of the electromagnetic transition amplitudes, and unrestricted Hartree- 
Fock method in atoms and nuclei violates rotational symmetry. Usually the symmetry is 
approximately restored in the final answer despite the violation on the way. We hope that 
the situation is similar in Takahashi's modified spin-wave theory. 

In our approach the dynamic structure factor Sm(w,<1) contains only two-quasiparticle 
intermediate states. It means that there is no simple vector excitation which would give 
CqS(u — f2 q ) contribution to the dynamic structure factor. In this sence the picture is 
similar to that in one dimensional Heisenberg chain where elementary excitation is spinon. 
We stress that the number of dynamic degrees of freedom is also similar to the S — 1/2 case 
because only doublets S z — ±1 are involved in the dynamics. However it would be probably 
wrong to say that the spin of elementary excitation S = 1/2, because still S z = ±1. More 
close analogy to the present situation is a gauge vector field in the axial gauge. In this sence 
the violation of the rotational symmetry is similar to the fixing of gauge. There are two 
separate questions in this situation 

1) Is the structure of spectrum in vector channel without narrow (^-functions valid for J\ — J 2 
model, or it is just a byproduct of the approximation? 

2) Is this technique applicable to t — J model? 

The answer to first question is not clear. I hope that further numerical simulations can 
elucidate this question. The answer to second question is more clear. In t — J model narrow 
structures in ^-distribution are much less important due to the very strong damping of spin- 
waves. In this situation the most important point is to avoid double counting of low-energy 
spin degrees of freedom. Takahashi's modified spin-wave theory definitely provides a correct 
count of these degrees of freedom. Therefore in my opinion this approach is reasonably 
justified for t — J model. 

I would like to thank M. Kuchiev, J. Oitmaa, R. R. P. Singh and A. Dotsenko for valuable 
discussions. This work forms part of a project supported by a grant of the Australian 
Research Council. 
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Figure caption, Fig.l 

Energy per site at J^j J\ = a ~ a c ~ 0.4. The value of a — a c is chosen in such a way that 
the gap at zero temperature equals u = 0.05. 
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Fig.1 



